Systematic study for particle transverse momentum asymmetry 
in minimum bias pp collisions at LHC energies 
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In PYTHIA6 (PYTHIAS) once the transverse momentum px of a generated particle is randomly 
sampled, p^ and Py are set on the circle with radius of pr randomly. This may largely suppress 
the development of the final hadronic state transverse momentum anisotropy from the initial state 
spatial asymmetry. We modify PYTIIIA6.4 by randomly setting p^ and py on the circumference 
of an ellipse with the half major and minor axes being pt(1 + 5p) and pt(1 — 5p), respectively. 
The modified PYTHIAS. 4 is then employed to systematically study the charged particle transverse 
momentum asymmetry in the minimum bias pp collisions at y^=0.9, 7, and 14 TeV. The ALICE 
data on the transverse sphericity as a function of charged multiplicity, {ST){Nch), are well reproduced 
with the modified PYTHIA6.4. It is found that the predicted charged particle V2 upper limit is a 
measurable value, ~ 0.12, in the minimum bias pp collisions at v^=7 TeV. We suggest a systematic 
measurement for the particle transverse momentum sphericity, eccentricity (ellipticity) , and elliptic 
flow parameter. 
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The event shapes are relevant to the properties of 
hadronic final state. One has employed the hadronic 
event shape measurements to test asymptotic freedom 
and to extract the strong coupling constant, etc. in the 
e~^e~ annihilation and lepton deep inelastic scattering, 
respectively, for a long time jl|, Q. Recently, the final 
hadronic event shapes in pp collisions at LHC energies 
have been measured by CMS ALICE 0,3, and AT- 
LAS 3, 13. 

In order to avoid bias from the boost along beam axis 
, this study is restricted to the transverse momentum 
plane. We start from the transverse momentum matrix 



of produced charged particles 



Pxi PxiPyi 
PviPx, 



Pi 



(1) 



where px^ is the transverse momentum of particle i, Px^ 
and Py^ are the transverse momentum components, and 
the sum here runs over the charged particles in a single 
event. The two eigenvalues of this transverse momentum 
matrix satisfy Ai + A2 = 1. If the two roots are ordered 
in Ai > A2, the transverse momentum sphericity is then 
defined as [51] 



St — 2A2. 



(2) 
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By intuitive construction one knows that St possesses 
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St 



(3) 



limits of 

0, pencil — like limit 

1, isotropic limit 
The event averaged transverse momentum sphericity is 
then denoted as (St), where (...) indicates an average 
over events. 

The dimensionless observable 

EM 



(4) 



(5) 



with limits of 

1 , pencil — like limit 
0, isotropic limit 
was first introduced in [9] to investigate the particle 
transverse momentum asymmetry. We are inspired by 
[lo| to identify this observable as transverse momen- 
tum eccentricity (ellipticity). According to the limits 
in Eqs. ^ and ([5]) one may intuitively assume that St 
is inversely proportional to ep, that is, a larger spheric- 
ity corresponds to a smaller ellipticity (eccentricity) and 
vice versa. Then the event averaged transverse momen- 
tum eccentricity (ellipticity) is 



Pv 



(6) 



\E^[Pl+Pl]/' 

Unfortunately, the transverse momentum eccentricity 
(ellipticity) was no longer a hot topic since then. 

The particle transverse momentum asymmetry was 
also proposed to be studied by the harmonic coefficients 
in the Fourier expansion of particle transverse momen- 
tum azimuthal distribution The second har- 
monic coefficient (elliptic flow parameter V2) is specially 
important because the large V2 of emitted particles is 
a characteristic feature of the hot and dense medium 
created in the ultra-relativistic nuclear collisions. It 
has contributed to the suggestion of a strongly coupled 
quark-gluon plasma (sQGP) observed in the nucleus- 
nucleus collisions at RHIC energies This elliptic 



flow parameter V2 is expressed 12| as 



{V2 



(cos[2(0 -*,)]), 



Px-Py 

,pI+pI 



(7) 



where 4> refers to the azimuthal angle of particle trans- 
verse momentum, stands for the reaction plane an- 
gle, and (...)p denotes an average over all particles in 
all events, which is a particle- wise average 17|. As men- 
tioned in [17], in order to take the multiplicity effect into 
account it is better to employ the event-wise average V2 
expressed as 



(cos[2(0 -*,)]), 



Px~Py 

,pI+pI 



(8) 



where (...)e indicates an average first over particles in a 
single event and then over events. Note that the event- 
wise averaged (w2)e is different from the (ep) in Eq. ([6|). 

As mentioned in ^ that the charged particle trans- 
verse momentum sphericity is measured (defined) "us- 
ing jets to represent the final state four-momentum." 
As the same elemental observables, and py. , appear 
in Eq. ^ and (|3]), the simultaneous measurement of 
sphericity and ellipticity (eccentricity) is not hard. This 
measurement is only effected by the jet reconstruction. 
However, the measurement of V2, whatever the event 
jlane method [12] or the Lee- Yang zero point method 
181 ] or the cumulant method is more model de- 



pendent, such as the decomposition of the nonflow 
etc. The cumulant method is further identified with 
two-, four-, and six-particle cumulants. The discrep- 
ancy among the measured V2 with different methods may 



reach 10-100% 



21 



22 



Recently, one even argued that 



231 . Therefore, we 



the event plane method is obsolete 
strongly suggest that the transverse momentum spheric- 
ity, ellipticity (eccentricity) , and the elliptic flow param- 
eter should be measured simultaneously, for the benefit 
of cross checking and reliable measurements. 

The ALICE data of charged particle transverse 
sphericity have been simulated with the different 
PYTIIIA6 tunes Q] and are best reproduced by the 
PERUGIA-2011 tune Q- In th^^ PERUGIA-2011 tune 
a number of theoretical parameters, relevant to the soft- 
and hard- processes as well as the parton distribution 
fmiction etc., are fitted to the data taken in hadronic 
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FIG. 1: (colour on line) Charged particle event averaged transverse sphericity as a function of charged multiplicity 
in the minimum bias pp collisions at (a) -^5=0. 9 TeV and (b) y/s—7 TeV. The ALICE data are taken from 5j and 

the PYTHIAS results are copied from [sf. 



decay at LEP, to the Tevatron min-bias data at differ- 
ent energies, to the Tevatron Drell-Yan data, and to the 
SPS min-bias data, etc. Therefore, it is possible that one 
parameter fitted to the hadronic Z° decay data, for in- 
stance, is good for pion production in nuclear collisions 
but not suitable for kaon production. The physics in 
PYTHIA6 for kaon production may even be biased. This 
situation, of course, may exist in other PYTHIA6 tunes 
quoted in Q|. This causes that the different PYTIIIA6 
tunes have different potentials describing a single ob- 
servable and that a single PYTHIAS tune has different 
potentials describing the different physics. It is also the 
reason that why the PERUGIA-2011 tune can be used 
with only PYTHIA6 but not with PYTHIAS [24[. Thus 



we prefer that the studied physics can be described by 
as simple modifications as possible to the existing model 
and/or by fitting as few model parameters as possible to 
the experimental data taken from the same reaction. 

In PYTHIA6 (PYTHIAS) once the transverse momen- 
tum pt of a generated particle is randomly sampled, 
and Py are then set on the circle with radius of px ran- 



domly. This strongly cancels the final hadronic state 
transverse momentum anisotropy developed from the ini- 
tial state spatial asymmetry. Therefore in this paper we 
modify PYTHIA6 (PYTHIA6.4 0) by randomly set- 
ting px and Py on the circumference of an ellipse with 
the half major and minor axes being J3t(1 + Sp) and 
Pt{1 — Sp), respectively. This change is also introduced in 
the particle/parton production process of hard scatter- 
ing, multiple interactions, initial- and final-state parton 
showers, as well as the adding of remnants 25|. Conse- 
quently, we re-calculate the pT — Px ^ Py ^^'^ modify 
the transverse momentum distribution of final hadronic 
states. It turned out that, provided 5p < 0.1 (a small 
perturbation) the change in the transverse momentum 
distribution is weak, so the comparison of transverse mo- 
mentum distribution with experimental data may not be 
effected. 

The modified PYTHIA6.4 is then employed to si- 
multaneously calculate the charged particle transverse 
sphericity, ellipticity (eccentricity), and elliptic flow pa- 
rameter in the minimum bias pp collisions at ^^=0.9, 
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7, and 14 TeV. In the calculations the model parameters 
are default 25| and the deformation parameter 6p intro- 



duced in this work is first assumed to be 0, 0.1, 0.2, and 
0.3, and then determined by fitting to the ALICE data on 
the event averaged charged particle transverse spheric- 
ity as a function of charged multiplicity, (S't)(A^c/i) fl- 
Note that the model reduces to the defauh PYTHIA6.4 
{25I when Sp=0. 

The calculated event averaged charged particle trans- 
verse sphericity as a function of charged multiplicity, 
{ST){Nch), is compared with the ALICE data [sj in 
Fig. [T] (a) and (b) for the minimum bias pp collisions at 
Y^=0.9 and 7 TeV, respectively. In the figure the solid 
circles are the ALICE data of "all" events taken from 
The red solid, green dashed, blue thick solid, and 
black dash-dotted curves are calculated by the modified 
PYTHIA6.4 with 5p^0, 0.1, 0.2, and 0.3, respectively. 
One sees in Fig. [T] (a) that the simulations reproduce 
the 0.9 TeV ALICE data with 6p=0 and 0.1. We also 
see in this panel that PYTHIAS is not so well compar- 
ing with the ALICE data, because it is not modified yet. 
Figured] (b) shows that the ALICE data of (S't)(A'^c/i) in 
the minimum bias pp collisions at TeV are better 

reproduced by the simulations with Sp=0.1. The cal- 
culated charged particle V2 upper limit is a measurable 
value of 0.119 (see Table U below). 

TableUgives the event averaged charged particle trans- 
verse sphericity, eccentricity (cUipticity), and the ellip- 
tic flow parameter in the minimum bias pp collisions at 
^/s=0.9 and 7 TeV. One sees in this table that: 

• The event averaged charged particle transverse 
sphericity decreases with increasing deformation 
parameter Sp but the event averaged charge par- 
ticle eccentricity (cUipticity) and the V2 upper 
limit (wj^ in the table) behave quite the oppo- 
site. Therefore, the above assumption that St is 
inversely proportional to ep (^2) is proved qualita- 
tively. We have to mention here that the V2 upper 
limit is calculated at the peak of the event-wise 



averaged V2[pt) hi |?;| <0.8. 

• The ALICE datum of event averaged charged par- 
ticle transverse sphericity in the minimum bias pp 
collision at -y/s=0.9 TeV is reproduced the best 
with (5p=0, consistent with Fig. [T] (a). Then the 
charged particle V2 upper limit is zero. 

• It appears that for the pp collision at y^=7 TeV 
the ALICE datum of event averaged charged par- 
ticle transverse sphericity, {ST)exp, can be repro- 
duced with Sp from to 0.2. 



As shown in Table HI the calculated (St) the is not so 
sensitive to Sp, unlike (ep) and V2^. {ST)the is also less 
sensitive to Sp than the charged particle event averaged 
transverse sphericity as a function of charged multiplic- 
ity {ST){Nch) in Fig.[T](b), since the former is obtained 
by averaging the later over Nch- However, considering 
Table H] together with Fig.[T](b) it is reasonable to choose 
(5p=0.1. So the charged particle V2 upper limit is a mea- 
surable value of 0.119, consistent with the existed predic- 
tion of 0.04 - 0.2 obtained by artificially setting param- 
eters in the models to upper limits for the pp collisions 



at Vs=7 and/or 14 TeV [26|-3- 

We then use Sp^O.l to predict the charged particle 
event- wise averaged W2('7) (0.5< pr <20 GeV/c) and 
V2{pt) (|?/| <0.8) in the minimum bias pp collisions at 
-ys=7 and 14 TeV as shown in Fig. [5] (a) and (b), re- 
spectively. We see in Fig. [5] that the charged particle 
V2 upper limit estimated at the peak of the V2{pt) in 
Fig. [2] (b) is larger than the one estimated at the peak 
of V2{rf) in Fig. [2] (a), because the former is counted in a 
rj interval of 1.6 but the later in an interval of 1. The V2 
upper limits, read at the peak of solid curve (for the pp 
collisions at y^=14 TeV) and the dashed curve (for the 
pp collisions at v^=7 TeV) in Fig. [2](b), are 0.194 and 
0.119, respectively. Thi s is con sistent with the existed 
predictions of 0.04 - 0.2 



26-2 



3- 



In summary, we have modified PYTHIA6.4 by set- 
ting the generated particles on the circumference of an 
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TABLE I: Event averaged charged particle transverse sphericity, eccentricity (ellipticity), and V2 upper Umit 
calculated at the peak of the event-wise averaged V2{pt) (Ivl <0.8) in the minimum bias pp collisions at y^=0.9 

TeV and y/s=7 TeV. 



Vs=0.9 TeV 



^=7 TeV 



''2 



Sp {St)Ixp {ST)thc (ep) «2 

0. 0.613±0.020 0.602 0.000 0.000 

0.1 0.597 0.109 0.105 

0.2 0.582 0.215 0.207 

0.3 0.555 0.315 0.302 



{ST)exp {ST)thc i^p) 

0.700±0.048 0.703 0.000 

0.696 0.109 

0.672 0.213 

0.632 0.310 



^2 

0.000 
0.119 
0.237 
0.347 



' obtained by averaging the ALICE data in Fig. [1] (a) over Nch and by quadratically 
summing the error bars in the ALICE data. 
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FIG. 2: (colour on line) (a) V2{ri) and (b) V2{pt) of charged particle calculated by modified PYTHIA6.4 with 
^p^O.l for the minimum bias pp collisions at ^/s=7 and 14 TeV, respectively. 



ellipse with the half major and minor axes being respec- 
tively Pt(1 + Sp) and pt(1 ~ ^p)i instead of on a circle 
with radius px originally. The modified PYTHIA6.4 is 
then employed to calculate the charged particle trans- 
verse sphericity, ellipticity (eccentricity), and the ellip- 
tic flow parameter in the minimum bias pp collisions at 
V^=0.9, 7, and 14 TeV systematically. By fitting to the 
ALICE transverse sphericity data a suitable deformation 
parameter 6p is obtained generally. With this Sp the el- 
liptic flow parameter as a function of 77 and px, V2{ri) and 
'V2{pt), are predicted for the minimum bias pp collisions 
at \/s—7 and 14 TeV. The V2 upper limits, estimated at 



the peak of the event- wise averaged V2{pt) calculated in 
\-q\ <0.8 with (5p=0.1, are 0.119 and 0.194 for the min- 
imum bias pp collisions at and 14 TeV, respec- 
tively. This theoretical prediction is consistent with the 



existed prediction of 0.04 - 0.2 |26|-3- 



We strongly suggest that the transverse momentum 
sphericity, ellipticity (eccentricity), and elliptic fiow pa- 
rameter are measured simultaneously, for the benefit of 
cross checking and reliable measurements. The theoreti- 
cal way employed in this work, that is, fitting 5p first to 
the measured charged particle transverse sphericity as a 
function of charged multiplicity and then employing the 
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fitted 6p to study the elliptic flow and higher harmonics, 
may be applicable to other works. 
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